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^^ Abstract: We propose the concept, synthesis, analysis and design of 

C*~) graphene-based plasmonic tunable low-pass filters operating in the THz 

band. The proposed structure is composed of a chemically-doped graphene 
^_^ strip deposit on a dielectric and a set of polysilicon DC gating pads 

located beneath it. This structure implements a stepped impedance low-pass 
filter for the propagating surface plasmons by adequately controlling the 
guiding properties of each strip section through graphene's field effect. A 



synthesis procedure is presented to design filters with desired specifications 
in terms of cut-off frequency and in-band and rejection characteristics. 

Ch The electromagnetic modeling of the structure is efficiently performed by 

combining an electrostatic scaling law to compute the guiding features of 
each strip section with a transmission line and transfer-matrix framework, 
which is also validated via full wave simulations. Excellent performance 
f^j is obtained for this frequency range in terms of roll-off, insertion losses 

£^ (rj 3 dB assuming a graphene relaxation time of T = 1 ps) and reconfig- 

uration capabilities (around 50%). These results, together with the high 
miniaturization associated with plasmonic propagation, are very promising 
for the future use and integration of the proposed filters with other graphene 
and silicon-based elements in innovative THz communication systems. 
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1. Introduction 

The unique electrical and optical properties of graphene and its ability to support ultra-confined 
low-loss plasmons are paving the road to the development of all-integrated graphene plas- 
monic devices at THz frequencies [1, 2, 3|. These devices may find application in different 
areas such as chemical and biological remote sensing, high resolution imaging and tomog- 
raphy, time-domain spectroscopy, atmospheric monitoring, and broadband picocellular or in- 
trasatellite communication networks [4, 5 6|. The desired integrated plasmonic structure must 
perform multiple complex functions, including signal detection, processing, and transmission. 
Currently, graphene-based plasmonic components ranging from waveguides [7 8 9 10|, anten- 
nas lTrnri"2 l. reflectoarrays [13], amplifiers fl4l . switches lfT31 . modulators [ 16|, reconfigurable 
frequency-selective surfaces ifTTl . or phase-shifters ifTSl to sensors [19, 20| have already been 
proposed and investigated. 

Terahertz systems also require filtering elements to select target frequency bands and re- 




Fig. 1 : Proposed graphene-based THz low pass filter of N' h degree. The structure consists of a 
monolayer graphene strip and N gating pads located beneath it. 



ject thermal radiation that may otherwise saturate sensitive detectors IF2T1I221 . Although sev- 
eral plasmonic guided filters have been developed at THz using noble metals and composite 
structures [23, 24l[25), graphene may lead to the development of novel filters with enhanced 
performance in terms of reconfiguration capabilities, miniaturization, insertion losses, and inte- 
gration with other silicon-based elements. However, the use of graphene for the design of such 
components has not been investigated yet. 

In this context, we propose the concept, analysis, and design of graphene-based plasmonic 
tunable low-pass filters. The structure is composed by a chemically-doped graphene strip and 
several independent polysilicon DC gating pads located beneath it. The strip supports the prop- 
agation of extremely-confined transverse-magnetic (TM) plasmons ||8] |9l whose guiding char- 
acteristics can be dynamically modified along the line by applying different DC bias voltages 
to the gating pads. In case that all pads are equally biased, the structure behaves as a simple 
plasmonic transmission line (TL) propagating the input waves towards the output port. When 
a different DC bias is applied to a gating pad, the guiding properties of the strip area located 
above are modified thanks to graphene's field effect. We apply this geometry to implement 
stepped impedance low pass filters, which are composed of a cascade of transmission lines al- 
ternating sections of high and low characteristic impedance. Importantly, the cutoff frequency 
of the filters can be dynamically tuned by simultaneously modifying the DC bias applied to the 
gates. A synthesis procedure is then presented to design filters with the desired cutoff frequency 
and in-band and rejection characteristics, directly providing the physical length of the gating 
pads and the required biasing voltages. The electromagnetic modeling of the structure is per- 
formed combining a transmission line model with a transfer-matrix approach. To this purpose, 
the recently introduced graphene electrostatic scaling law |9] is applied to efficiently compute 
the propagation constant of the modes supported by the strip as a function of its width, sur- 
rounding media, and applied electrostatic DC bias. The proposed approach allows the accurate 
analysis of the desired filters in just seconds, avoiding large simulation times of general purpose 
full-wave software. 

In order to illustrate these concepts, several low-pass filters are designed and analyzed, eval- 
uating their performance and reconfiguration capabilities in the THz band. These initial results 
are very promising for the integration of graphene-based filters in future THz plasmonic devices 
and communication systems. 
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Fig. 2: Equivalent transmission line model of the graphene-based filter shown in Fig 111 

2. Proposed structure 

The proposed structure, depicted in Fig. [T] comprises a chemically doped graphene ribbon 
deposited on a dielectric substrate, and a number of polysilicon gating pads beneath the strip. 
Graphene can be characterized by a complex-valued surface conductivity a, which can be mod- 
eled using the well-known Kubo formalism [26], and mainly depends on graphene chemical po- 
tential fi c , relaxation time T and temperature. One of the most interesting features of graphene 
is that its chemical potential can be tuned over a wide range (typically from -1 eV to 1 eV) by 
applying a transverse electric field via a DC biased structure, such as the one proposed here. An 
approximate closed-form expression to relate ji c and the applied DC voltage (Voc) is given by 



^J 7 ^, (1) 

where C ox = £ r £o/t is tne S ate capacitance, and e r and t are the permittivity and thickness 
of the gate dielectric. Promising measurements have recently been reported in E71l . where a 
variation in the chemical potential from 0.2 to 0.6 eV was obtained by applying a bias voltage 
from 0.2 to 1.73 V. 

In addition, graphene monolayers support the propagation of surface plasmons polaritons 
(SPP) at THz frequencies with moderate losses and extreme confinement. Several authors have 
studied the characteristics of the SPPs propagating along graphene ribbons OH |9), and trans- 
mission line models have been successfully utilized to describe this type of structure irT4ll28l . 
Using this approach, the structure of Fig.fTlcan be modeled as a cascade of transmission lines, 
as shown in Fig. [2] The complex-valued characteristic impedance and propagation constant of 
the transmission lines depend on graphene's chemical potential and will be largely modified by 
the DC voltage applied to the gating pads, allowing the synthesis of the proposed filters. 

3. Synthesis and Modeling of Graphene-based Filters 

3.1. Synthesis Procedure 

The goal of this section is to design a lowpass filter using the structure of Fig.fT] with the desired 
cutoff frequency, rejection characteristics and inband performance. This structure implements 
a so-called stepped impedance lowpass filter (29), whose equivalent network is presented in 
Fig. [2] As seen, it is composed of the connection of N transmission line sections with lengths 
Ik, propagation constants jk and characteristic impedances Z^. 

The design procedure starts with the calculation of transfer and reflection polynomials to 
satisfy the desired specifications in terms of in-band and rejection characteristics (order of the 
filter, rejection level, return losses) 

*»<»>= i*W. Sn=F M (2) 
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Fig. 3: Characteristic impedance (a) and normalized phase constant (b) of the first waveguide- 
like surface plasmon propagating along different graphene strip configurations versus the chem- 
ical potential jj. c . Solid lines have been obtained using the electrostatic approach described in 
Section III and markers indicate values computed with HFSS. Graphene parameters are T = 1 ps 
and J = 300 K. 



The calculation of the P(s) and E(s) polynomials is done analytically for most useful trans- 
fer functions, including Butterworth and generalized Chebyshev responses. Some useful tech- 
niques are reported in ||29l . The next step in the design procedure is the election of the electrical 
length 9 C of the individual line sections. This parameter directly determines the periodicity of 
the frequency response when the ideal polynomials are implemented with transmission lines. 
Smaller values of 9 C result in a wider spurious-free range, while requiring more abrupt changes 
in the line impedances. Having decided the value of 9 C , a recursive technique [29 1 is applied 
to extract the normalized values of characteristic impedances (Z^). Then, the de-normalization 
of the calculated characteristic impedances to the real port impedances (Zo) used in the filter 
implementation is done as Z# = Z^Zq. 

The final and crucial design step consists in finding the design parameters of the physical 
structure in Fig. IT] to implement the prototype circuit at the desired cut-off frequency. To this 
purpose, the required de-normalized impedance values (Zj) obtained during the above proce- 
dure are synthesized using the SPP properties of graphene strips. This can be most efficiently 
accomplished by adjusting appropriate chemical potentials (ji c ) through electrostatic gating. As 
an illustrative example of this electrical control, we show in Fig. [3]the real part of the charac- 
teristic impedance and normalized phase constant of the first waveguide-like mode propagating 
along a graphene strip, computed versus the chemical potential for graphene waveguides of dif- 
ferent characteristics. In our work we have decided to maintain fixed the width of all the strips, 
and to adjust the impedances of all sections (Z^) by acting only on the chemical potential by 
applying an electrostatic gating. It can be observed in the figure that, for several strip widths and 
frequencies, the impedances that can be obtained vary from around 20 k£2 to several hundreds 
of ohms for practical chemical potentials in the range 0—1 eV. Note that the calculation of 
the chemical potentials for all sections, following this procedure, also fixes the values of their 
propagation constants (j3^). This information, and the election of {9 C ) done during the synthesis 
phase, allows the calculation of the physical lengths of all graphene strip sections in Fig.[T]with 



the straightforward relation 

k ~m (3) 

Note that this synthesis procedure is strictly valid when lossless transmission lines are con- 
sidered. The presence of losses in the real structure leads to small deviations between the actual 
filter response and the expected synthesized function l29l , with some extra round-off inside the 
filter passband. 

3.2. Electromagnetic Modeling 

The electromagnetic modeling of the structure shown in Fig. [1] is based on the analysis of its 
equivalent network (see Fig. [2} using a transmission line formalism combined with an ABCD 
transfer-matrix approach [30]. The main difficulty of this approach is to accurately compute 
the complex propagation constant and characteristic impedance of the surface plasmons propa- 
gating along each section of the strip. The properties and dispersion relation of these plasmons 
have been studied by several authors EH). Unfortunately, the desired dispersion relation does 
not admit an analytical solution and one has to resort to full-wave numerical techniques. Con- 
sequently, the use of the transmission line approach requires the numerical analysis of several 
graphene strips, which is a tedious and time-consuming task and does not incorporate any ad- 
vantage over performing the full-wave analysis of the whole structure. 

In this context, |9| studied the electrostatic nature of surface plasmon propagating along 
graphene strips and proposed a simple scaling law able to obtain general dispersion relations. 
This electrostatic approach, which assumes that the strip width is much smaller than the wave- 
length, establishes that plasmon properties are solely determined by the strip width (W), sur- 
rounding media (e r ), and graphene conductivity (<r). Then, once the propagating features of 
a given plasmonic mode have been obtained, they can be scaled to any arbitrary situation (in 
terms of W, £,- and a) by using the scaling parameter 

\m\a(fa)\ 

where % = 2/(e, + 1) is a constant that models the dielectric media and fp is the frequency 
where the surface plasmon propagates with a phase constant /3 . Note that, due to the electro- 
static approach employed to derive Eq. |4](see |9|), the scaling parameter 77 is independent of 
the operation frequency. 

This scaling law is applied to efficiently compute the SPP propagating features along any 
graphene strip. The process is as follows. First, the scaling parameter 77 related to the desired 
surface mode is obtained. To this purpose, a single run of full-wave simulations is required 
to obtain the dispersion relation of the mode, thus allowing to compute rj using Eq. [4] Fig. [4] 
confirms that rj only depends on the product fiW [9| and that surface plasmons propagating 
on different strips leads to the same scaling parameter. Second, the phase constants of surface 
plasmons propagating along strips of arbitrary characteristics are computed using tj. To this 
purpose, the scaling parameter is computed at the desired operation frequency fa and strip 
width W using Eq.Hand then the corresponding value of fiW is interpolated using the informa- 
tion of Fig. [4] Since the strip width W is known, the phase constant of the surface plasmon is 
easily obtained. Finally, the attenuation constant is found as a = l/(2L p ), where L p is the 1/e 
decay distance of the power. The propagation distance of plasmons is mainly controlled by the 
electron relaxation time of graphene T, and it can be approximately obtained by v g T, where v g 
is the plasmon's group velocity (v g = dOp/dj3) Q. Therefore, the attenuation constant may be 
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Fig. 4: Scaling parameter r\ versus the product of the phase constant and the strip width (j5W) 
for different graphene strip configurations. Graphene parameters are T = 1 ps and T = 300° K. 
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Once the complex propagation constant 7=0! 
its characteristic impedance is obtained as QUI 



jj3 of the propagating TM plasmon is known, 
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where E e ff = (1 + £ r )/2. The accuracy of this approach is validated in Fig. [3] where the char- 
acteristic impedance and phase constant of plasmons propagating along different strip configu- 
rations are computed using both the proposed technique and the finite element method (FEM) 
software HFSS. There, the graphene strip is modeled as an infinitesimally thin layer where 
surface impedance boundary conditions are imposed (Z mr f = 1/(7.) 

The combination of the graphene strip scaling law with a transmission line and transfer- 
matrix approach permits an extremely fast electromagnetic analysis of the proposed filtering 
structure (see Fig.[T|, allowing an efficient implementation of the synthesis technique described 
in the previous section. For validation purposes, the electromagnetic response of the complete 
proposed structure is also obtained using the full-wave software package HFSS. 

4. Results and discussion 

In this section we design and analyze two low-pass filters which are implemented using the 
structure shown Fig.fTI For the sake of generality, the filters have been designed to have different 
order and cutoff frequencies, considering various strip widths and dielectrics. The performance 
of the filters is presented in terms of their scattering parameters, referred to the characteristic 
impedance of the graphene sections at the input and output ports. A comparison between the 
transmission line model combined with the scaling law and full-wave results using HFSS is 
shown for both cases, validating the accuracy of the proposed electromagnetic modeling. In 
addition, the reconfiguration capabilities of the filters are investigated in detail. In our study, we 
consider graphene with a relaxation time T of 1 ps and a temperature of T — 300° K. 



Table 1: Design parameters of the first example: a 7th degree filter. 
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Fig. 5: Scattering parameters of a 7 r/ ' degree filter implemented using the structure depicted 
in Fig. [T] The design parameters are shown in Table [T] the strip width is 150 nm and ion gel 
of e, = 1.8 is employed as a dielectric, (a) Nominal filter designed to have a cutoff frequency 
of 2.3 THz. Results are obtained using the ideal synthesis procedure, the transmission line ap- 
proach combined with the scaling law, and the commercial software HFSS. (b) Reconfiguration 
possibilities of the filter obtained by adequately controlling the DC bias of the different gating 
pads. 



In the first example, we consider a graphene strip of 150 nm deposited on an ion gel dielectric 
of permittivity 1.8 Il27l . Using this host waveguide, we have designed a 1 th degree filter with 
a cutoff frequency of 2.3 THz. The Chebyshev polynomials were computed using standard 
techniques |29| for a maximum theoretical return loss of 30 dB, and the electrical length 9 C 
at the cutoff frequency was set to 37°. These values were chosen to yield practical values of 
chemical potential through the synthesis procedure explained in the previous section. The final 
design parameters of the filter are shown in Table [1] Note that filters of this degree with better 
roll-off characteristics could be synthesized, but this would further increase the required range 
of chemical potential values achievable in the structure. 



Fig. 5a shows the frequency response of the filter, computed using the proposed transmission 
line approach and the full-wave commercial software HFSS. A good degree of agreement is 
observed, with very similar attenuation profile and average level of return loss in the passband. 
Moreover, the structure presents a low insertion loss level, around 3 dB, which is a remarkable 
value at this frequency range. Insertion losses increase for lower values of graphene's relaxation 



Table 2: Design parameters of the second example: a 9th degree filter. 



Section 


Z 


/(nm) 


t^c, nominal \P * / 


Vc.tunedl (eV) 


detuned! (eV) 


Ports 


1 


200 


0.17 


0.35 


0.46 


1,9 


1.36 


156 


0.1 


0.21 


0.274 


2,8 


0.58 


367 


0.43 


0.93 


1 


3,7 


2.24 


94 


0.035 


0.01 


0.12 
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0.44 
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0.69 
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Fig. 6: Scattering parameters of a 9 th degree filter implemented using the structure depicted 
in Fig. [T] The design parameters are shown in Table 13] the strip width is 100 nm and quartz 
of e, = 3.9 is employed as a dielectric, (a) Nominal filter designed to have a cutoff frequency 
of 3.3 THz. Results are obtained using the ideal synthesis procedure, the transmission line ap- 
proach combined with the scaling law, and the commercial software HFSS. (b) Reconfiguration 
possibilities of the filter obtained by adequately controlling the DC bias of the different gating 
pads. 



time, but due to the short length of the structure in terms of the plasmon wavelength, they 
remain at acceptable levels (w 5 dB for a relaxation time of T = 0.5 ps). The slight differ- 
ence in the cutoff frequency between both approaches is due to the mono-modal nature of the 
transmission line approach, which neglects the possible excitation of higher order modes at the 
connection between two adjacent sections. This phenomenon modifies the phase condition of 
the circuit, which is now fulfilled at a slightly different cutoff frequency. Interestingly, this ef- 
fect appears to be uniform along the structure, allowing to adjust the filter's cutoff frequency 
by a small overall scaling of the section lengths. Fig.|5b]illustrates the reconfiguration possibil- 
ities of the designed filter. By adequately controlling the DC voltage applied to the gating pads, 
following the synthesis procedure described in Section III, the overall electrical length of the 
device can be modified thus leading to an electric control of the filter's cutoff frequency. 



For the sake of clarity, only two different possible reconfiguration states are shown in the 
figure (their corresponding design parameters are shown in Table [TJ. However, a continuous 
range of cutoff frequencies can be easily synthesized. Importantly, the proposed filter allows a 
dynamic control of the cutoff frequency over 50%. This value is mainly limited by two factors. 
First, the fixed physical length of the gating pads imposes a limit to the maximum frequency 
shift attainable while maintaining an acceptable variation to the attenuation profile and return 
loss. And second, the values of the chemical potential employed are limited, by technological 
reasons ETl . to the range of — 1 eV. 

The second example is based on a waveguide of different characteristics, composed of a 
strip width of 100 nm deposited on a quartz substrate of £, = 3.9. The filter's degree has been 
increased to N = 9, and the cutoff frequency is set to 3.3 THz, with 9 C = 39°. The complete 
design parameters of the filter are shown in Table [2] and the frequency response is plotted in 
Fig. [6a] Very similar conclusions compared to the previous design can be drawn from the filter 
response. Importantly, the difference in the cutoff frequency between the two electromagnetic 
analysis is also of similar relative magnitude to that of the previous example. Finally, the tunable 



features of this filter are shown in Fig. 6b where around 50% of cutoff frequency dynamic 
control is obtained. 

5. Conclusion 

We have proposed and analyzed graphene-based plasmonic low-pass filters in the THz band. 
An efficient synthesis method has been presented, based on classical microwave filter theory 
adapted to graphene SPPs through a recently introduced electrostatic scaling law. Several re- 
sults, verified with full-wave simulations, showcase the filtering capabilities of the proposed 
structure and its enhanced performance compared to the alternatives in this frequency range, 
in terms of roll-off, miniaturization and real-time reconfigurability. In particular, through ade- 
quate use of electrostatic gating, the cut-off frequency of the device can be continuously tuned 
over a wide frequency range, resulting in a degree of tunability not achievable with any other 
existing technology in the THz band. We envision that the proposed reconfigurable filtering 
structures may be useful in future all-integrated graphene plasmonic THz devices, sensors and 
communication systems. 
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